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SC 29208 

A systematic formalism for quantum electrodynamics in a classical uniform magnetic field 
is discussed. The first order radiative correction to the ground state energy of an electron 
is calculated. This then leads to the anomalous magnetic moment of an electron without 
divergent integrals. Thorough analyses of this problem are given for the weak magnetic field 
limit. A new expression for the radiative correction to the ground state energy is obtained. 
This contains only one integral with an additional summation with respect to each Landau 
level. The importance of this formalism is also addressed in order to deal with quantum 
electrodynamics in an intense external field. 

§1. Introduction 

The discovery of the Dirac equation was essential to explain atomic phenomena. 
However, it turned out that nature requires modifications of the Dirac equation 
in several situations. The first difficulty was in explaining the phenomenon called 
“Lamb shift.” According to the Dirac equation, the energy levels 2 Si/ 2 and 2Pi/ 2 
of the hydrogen atom are degenerate. However, experiment shows that these states 
are separated by 1057Mhz. This experimental fact cannot be explained purely from 
the Dirac equation in a Coulomb field. 1 ) There were many attempts to resolve this 
difficulty based on phenomenological arguments, but none of them could explain 
from the first principle. It was only after the establishment of quantum field theory 
that a satisfactory explanation to this problem was given. Thus the quantum field 
theory of interactions between electrons and photons, quantum electrodynamics, was 
born and has been followed by great sequent successes. 

Quantum electrodynamics is one of the most successful idealized theories in 
physics. Agreements between theories and experiments are getting closer and closer 
with great precision. For example, the gyromagnetic factor of the electron, which 
is 2 according to the Dirac equation, is evaluated by one Feynman diagram in 1st 
order perturbation, 5 diagrams in 2nd order, 72 diagrams in 3rd order, and 891 di¬ 
agrams in 4th order. Every higher order correction result justifies the precision test 
for quantum electrodynamics. There is no doubt this theory is correct. However all 
these results and all other computations higher than the second order perturbation 
involve the following mathematical difficulties. If one computes physical quantities, 
such as energy of the electron, the polarization of photons, etc., one encounters inte¬ 
grals which diverge typically in the ultraviolet region. Physicists are more optimistic. 
We solved this infinity problem with the help of renormalization techniques. Renor- 
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malization can extract finite results from infinities, which agree with experimental 
numbers. Although there is no logical reason and justification to adopt this working 
rule, quantum field theory cannot be applied without renormalization. Moreover, we 
need the remornralization procedure even in low energy physics, such as condensed 
matter physics. 

We should be satisfied with the present success of quantum field theory to some 
extent, and go further to solve more problems in natural phenomena. We should 
however, sometimes look back and ask ourselves fundamental questions. 

Dirac was probably the most well known physicist with pure soul. He insisted 
that this renormalization be considered as a temporary developing stage of theo¬ 
retical physics. His attitude toward this problem can be seen from his works, for 
instance “Quantum Electrodynamics without Dead Wood”, 2 ) and his last article 
“The Inadequacies of Quantum Field Theory”. 3 ) Unfortunately, Dirac could not 
solve this problem completely. 

In this paper we would like to revisit this problem from a different point of view. 
One of the key points in this problem is to recognize that we are dealing with infi¬ 
nite dimensional Hilbert space in quantum field theory. Associated with an infinite 
degrees of freedom, representations of the canonical (anti-)commutation relations 
cannot be uniquely determined. It was however, rarely used to explain physical 
phenomena until the work of Umezawa. 4 )’ 5 ) His simple but elegant observation con¬ 
cluded that each different representation corresponds to different phases of physical 
systems. In a superconductor for example, superconducting states occur by a tran¬ 
sition from one representation to another representation accompanying spontaneous 
breakdown of the gauge symmetry. Here one representation corresponds to normal 
conducting states and the other superconducting states. Thus the same Hamilto¬ 
nian describing electrons in matter can describe different physical states. Therefore 
it is crucial to choose the proper representation of the canonical (anti-)commutation 
relations concerned, or in other words, the proper choice of the Fock space. Since 
the Fock space is constructed by a cyclic operation of creation operators on the Fock 
vacuum, the correct Fock vacuum has to be chosen from infinite many vacua. 

From Umezawa’s point of view, renormalization is understood as a procedure of 
changing the representation of operators. Since we have to stick to one representa¬ 
tion throughout the calculations, a wrong initial choice of Fock space may result in a 
renormalization procedure. To be specific, if one starts from a non-interacting parti¬ 
cle Fock space, it is necessary to carry out the renormalization procedure in order to 
calculate physical quantities for interacting particle systems. This observation seems 
to imply that any perturbation in quantum field theories require renormalization. 

With these observations it is quiet natural to apply the above idea to quantum 
electrodynamics. We will particularly consider the case where an electron is inter¬ 
acting with photons within a classical uniform magnetic field. The Lagrangian of 
the system is written as 

k-'tot, — L e i T Lem T -£mt, (1) 

where L e i, Lem, and L rn t are the Lagrangian for bare electrons, bare photons, 
and the interaction between electrons and photons respectively. In the presence of 
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classical fields, the photon term is written 

C-em = C^em + & 


cl 

- 'EM j 


( 2 ) 


where F‘e M represents the Lagrangian for the given c-number classical fields. Hence, 


r _ r , I r ( l rn _i_ rcl . nqm . r 

E-tot ~ E e l T E fm -f- Eem + E- f Z- 


( 3 ) 


The basic idea is to treat the first term and the last term together and quantize the 
electron field in a background of classical fields. Namely, we will consider 

C tot = l el + Cr M + CZ, ( 4 ) 

where C e i is the Lagrangian of electrons in classical fields, and E l '"' t is the interaction 
between electrons in classical fields and photons. 6 ) 

A computation of the anomalous magnetic moment of the electron using this 
idea was, to our knowledge, first done by Luttinger in 1949. 8 ) His observation was 
different from us, but he obtained the same first order perturbation result as is found 
with other methods but without divergent integrals. It seems that his work got little 
attention at the time. The main reason seems that his calculation was not relativis- 
tically covariant. Therefore, it was difficult to generalized to other processes. Later 
in his formulation of the proper time method, Schwinger independently obtained the 
same result. 9 ) Based on Schiwinger’s proper time method, Schwinger’s collaborators, 
and Ritus and his collaborators extensively studied the above idea, i.e. quantizing 
the Dirac field in classical background fields. 10 ) Their main interest is in applications 
to quantum electrodynamics in intense external fields, which cannot be handled in 
ordinary quantum electrodynamics. 

The main objective in this paper is to confirm Luttinger’s work from the mod¬ 
ern language, i.e. diagram techniques, and to establish the basic formulation in this 
problem. We would also like to discuss the relation between his result and others’ 
one. By reexamining this problem, we give an alternative expression for the first 
order radiative correction to the ground state energy. The final expression contains 
only one integral and a summation over the Landau levels. This new expression 
may be more useful than the previous result particularly for numerical evaluation. 
The physical interpretation is also clear by expressing each Landau level contribu¬ 
tion separately. Another motivation to revisit this classical problem is from recent 
astrophysical observations of extremely strong magnetic field, such as magnetars and 
pulsars. The strength of magnetic fields in these objects are much above the critical 
field in quantum electrodynamics, i.e. B c = m 2 c 3 /(\e\h) — 4.4 x 10 13 G. As will 
be discussed in this paper, the formalism present is essential in order to deal with 
strong external magnetic field problems. 

The organization of this paper is as follows. In order to make our paper self- 
contained, first few sections are devoted to give a review on quantization of the 
Dirac field with classical c-fields. The readers who are familiar with this can skip 
sections 2~4 without any problems. Exact solution of the Dirac equation in a uniform 
magnetic field is reviewed in Sec. 2. Notations and conventions are also summarized 
in Sec. 2. The Dirac field is then quantized using this complete orthonormal basis 
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in Sec. 3. The modified Dirac equation due to the radiative correction is discussed 
based on the mass operator formalism in Sec. 4. One loop correction to the mass 
operator is evaluated in Sec. 5 and Sec. 6 . The anomalous magnetic moment of an 
electron is calculated in a one loop level in Sec. 7. Divergent integrals appeared in 
this formalism is discussed in Sec. 8 . In Sec. 9 a compact expression is given for the 
energy shift within the lowest Landau level approximation. The paper closes with 
conclusion and discussion in Sec. 10. Throughout the paper we set h = 1 and c = 1. 

§2. Exact Solution of the Dirac equation 

The exact solution of the Dirac equation in a uniform magnetic field is well 
known. 13 ) To establish notations we give a brief review on the exact solution based 
on the operator formalism. The Dirac equation in an electromagnetic gauge field A^ 
is 

(ijl-m) |V>) = 0, (5) 

where Jjl = 7 ^77^ t and 77 M = P M — eA^. The Dirac gamma matrices 7 ^ satisfy the 
anti-commutation relation : 

{ 7 m ! Ivj = lv.lv + Ivin = 2 g^I, ( 6 ) 

where g jiv = diag(l — 1 — 1 — 1) is the Minkowski metric, and 7 is the 4x4 unit 
matrix. The electromagnetic field tensor F^ u is defined by the commutation relation 
between 77 M as 

[77^,77^] = -ieF^ w . (7) 

Here a charge of a particle is denoted by e including its sign, for instance e = — |e| < 0 
for the electron. In eq. (7) the canonical commutation relation [P^,X U \ = ig^ u is 
used, and the gauge field A M is assumed as a function of X (In this paper we 
express the position and momentum operators as X^ and P M respectively. Their 
eigenvalues are denoted as x M = ( t,x,y,z ) and = (j>t,Px,Py,Pz)-) I n the coordi¬ 
nate representation, the momentum operator is expressed as P^ = id/dx^. Then, 
the electromagnetic field tensor F^ u reads a familiar form : 

F^u = d^Av - dvAp. ( 8 ) 

By definition F^ is anti-symmetric, i.e. F^ = —F Ufl . 

Consider a uniform magnetic field where 77* = (0,0,77) along the z-axis with 77 
a constant. The non-vanishing components of the field tensor F^ u are then 

F‘2\ = —F\2 = 77, (9) 

and all other components are zero. To solve the first order (in P /t ) Dirac equation 
in external fields, it is convenient to bring it into the second order equation from 
multiplying by Jjl -\-m from the left, 

_ m 2_ 


{$ + m)(l/l - m) |V>) = (77 M 77 M 


( 10 ) 
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Here o^ u = ?'[7^,7 I/ ]/2 and eq. (7) is used. For the uniform magnetic field, the last 
term in eq. (10) is 

= - eB ( o ■ < u > 

Therefore it is possible to decompose the Dirac four spinor into each component 
which is also an eigenspinor of the z-conrponent of the spin matrix F 3 . Define the 
orthonormal spinors w a (a = 1,2, 3,4) by 









/°\ 


0 


1 


0 


0 

Wl = 

0 

, W2 = 

0 

, m = 

1 

, W A = 

0 


w 


\V 


W 


vJ 


These four spinors are eigenspinors of the spin matrix F 3 : F 3 w a = o a w a where 
eigenvalues o a are 

^ = +1(0 = 1,3), o a = -1 (a = 2,4). (13) 

Expanding the four spinor by the eigenspinors w a as | if) = w a\fa), the second 
order equation for the four spinor is now reduced to the one for scalar functions | f a ) 

(n^n^ - m 2 - o aQ 2 /2)\<j> a ) = 0. (14) 

where g = \/—2eB with the assumption — eB > 0, i.e. the direction of the constant 
magnetic held is upward along the z-axis. Since the only non-vanishing commuta¬ 
tion relations among 17^ is [TTi, 7 T 2 ] = ieB = —ig 2 /2, it is possible to reduce the 
problem to a two-dimensional one by separating these degrees of freedom. Choose 
the symmetric gauge for the gauge held : 

= (0, -BX 2 / 2, BX 72,0). (15) 

In this gauge IIq = Pq and IT 3 = P 3 . Then, the eigenfucntions for these operators are 
just plane waves with eigenvalues pt = E and p z respectively, e.g. P 3 \p z ) = Pz\Pz)- 
Therefore, essentially the problem is to solve the following equation on the xy-plane 

(n 2 + n 2 +m 2 + o a6 2 /2+p 2 z - E 2 )\(t>' a ) = 0 , (16) 

where \</> a ) = \E)\cf)' a )\p z ). 

For later convenience, we solve eq. (16) with the aid of the complex variable 
representation (holomophic representation) for ITi and I 72 - 14 ^ Introduce a complex 
variable £ = g(x + iy)/{ 2 \/ 2 ) and its complex conjugate f = g(x — iy)/{2\/2) in the 
coordinate representation. The differential operators with respect to £ and £ are 
related to d x = —iP\ and d y = — iP -2 by 
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The covariant differential operators D± = —ill\ and D 2 = —HI 2 are then expressed 
in terms of these new complex variables £, d, and d : 


D -^ (8 + 5 ) - 571 K - f) - 
D2 = i ^. (B ^ B)+ '^ {(+ ^ 


Hence, 

n\ + n 2 = -^ 2 (233 - 2 &-d£ -£d + B^ + ^B). 

Finally, introduce two sets of creation and annihilation operators by 

a t = (-3 + 0/V2, a = (3 + 0/V2, 
ft = (-3 +0/V2, (3=(3 + Z)/V 2. 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 


They satisfy the canonical commutation relations : 

[a,a l ] = [(3,ft] = 1, (24) 

and all other commutations are equal to zero. It is easy to show 

Uf + n% = Q 2 (a*a + -), (25) 

and this does not depend on the other set of creation and annihilation operators fv 
and (3. Then eq. (16) can be solved algebraically by 

[g 2 cYa + m 2 + p 2 + (cr a + 1 )g 2 /2 - E 2 ]\<B' a ) = 0. (26) 

Following an elementary exercise, the eigenvalues of ata are labeled by integers. 
Since a a = ±1, the energy eigenvalues are also labeled by integers n : 

E n = ±y 'm? + p 2 + ng 2 = ±e n (n = 0,1,2, • • •). (27) 

The ground state |0) of the system is defined by 

or|0) = /3|0> = 0. (28) 


The Hilbert space is then constructed by a cyclic operation of at and /?t on |0), viz., 


(a^) n (0^Y 

|n,l)= l >4 j |0) ( n , ^ = 0,1, 2, • • •). (29) 

Vn'.il 

These states consist of a complete orthonomal basis : 

(n,£\n ,£') = 5 n ySt/'. (30) 

The completeness relation is satisfied, 


^2 \n,e)(n,e\ = i, 

n,i =0 


(31) 
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where I is the identity operator. Using this basis, the energy eigenstates correspond¬ 
ing to the energy E n are found as 

W a { n ^)) = \ n ~ M> ( a = !,3), (32) 

|</4(M)) = \ n A) (a = 2,4), (33) 

where the ground state (n = 0) is realized only for a = 2,4. To treat the ground 
state on the same footage as the other excited states, it is convenient to use the 
convention : 

| - M) = 0 for all t. (34) 

Therefore the stationary solution of second order equation (I7 M I7 M — m 2 —a a g 2 /2)\(f a ) = 
0 is 

\(f> a (n,£,p z ;t)) =e~ tEnt \n-l,£)\p z ) (a = 1,3), (35) 

| <p a (n,£,p z ;t)) = e~ lEnt \n,Z)\p z ) (a = 2,4). (36) 

A similar convention will be used : 


\<j) a (—l,£)) = 0 for all £. 


(37) 


The solution of the original first order Dirac equation can be found as follows. 
Let | if) be the solution of second order Dirac equation : 

(ijl + — m) |y!) = 0, (38) 

then \if) = (Jjl + m)\ip) is the solution of the first order Dirac equation : 

(]/[ — m)\if) = (ijl — + m) |V>) = 0. (39) 

This argument becomes empty when \if) is also a solution of the first order equation, 
i.e. — m)\il>) = 0. It is straightforward to obtain an explicit representation of the 
operator ifl + m as follows : 


and hence, 


^IJi + 7 2 TI-2 = i^[a(7 1 + ry 2 ) - a t (7 1 - Z7 2 )] 



/ 0 0 

0 

a\ 



0 0 

— oft 

0 


= ig 

0 —a 

0 

0 

1 


0 

0 



fm + P 0 

0 

Pi 


iga \ 

0 

m + Pq 

—iga' 


-Pi 

-P-3 

—iga 

m — P 0 


0 

\ iga t 

Ps 

0 

m - P 0 J 


(40) 

(41) 


(42) 


The following standard representation of the gamma matrices are used in the above 
expressions : 


7 ° = 


7 = 


I 0 
0 -I 


0 a 1 
-a 1 0 


(43) 
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where I is the 2x2 unit matrix and a 1 are the Pauli matrices. The positive energy 
solution is obtained by applying the operator l/l + m to w\\4>i) and with 

Sin — d~£rj> 


\i>i{n,l,p z )t)) = Ne l£nt \p 


((e n + m)\n - 1,£)\ 

0 

p z \n - l,£) 

\ iy/ng\n, £) J 

( 0 \ 

(e n + m)\n,£) 

—iy/ng\n — 1,7) 

\ —p z \n — 1,7} ) 

The negative energy solution is obtained from and with E n = 


\^ 2 {n,£,p z ;t)) = Ne l£nt \p. 


(44) 


(45) 




\h{n,t,p z \t)) = Ne +l£nt \p z ) 

1^4 (n,£,p z ;t)) = Ne +l£nt \p z ) 


( p z \n — 1, 7) \ 

i\/ng\n, £) 

(e n + m)\n— 1,7) ’ 

V 0 / 

(—i\fng\n - 1,7)\ 
-p z \n - 1,7) 

0 


\ (£n + m)\n,£) J 


(46) 


(47) 


Except the ground state, each energy state is doubly degenerate with respect to the 
spin states of second order equation. Each energy level is also infinitely degenerate 
with respect to integers £. The normalization constant N is determined by 

('ipa{n,£,Pz)\^a(n,£,p z )) = 2N 2 £ n (e n + mi) = 1 (for each a). (48) 


Hence, 

N = [2e n (e n + m)] -1/2 . (49) 

This solution is a complete orthonormal, 

{i’a(n,£,p z )\^a(n',£',p' z )) = 6 a ,a'S n ,n'6e,p6(pz ~ p' z ), (50) 

and the completeness relation holds : 

^2 Jdp z \ipa{n,£-,t,pz))(ipi(n,£-t,p z )\ = /, (51) 


by use of relations (31) and f dp z \p z ){p z \ = I- 

Note that explicit representations for the x, y variables have not been obtained 
so far. The explicit xy-coordinate representations are constructed by a projection 
of (x, y\ onto the above solution. Equally, an ^-representation can be obtained by 
a projection of (£,£|. In the latter case, the ground state /(£,£) = (£, £|0) is easily 
found from two simple differential equations : 

(d+ <£)/(£,£) = 0, (B +Of (£,£) = 0. 


(52) 
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The solution is 

/(£,£) = y^-ex :p(-££). (53) 

Here the normalization constant is determined by 

Jdii(o m,£)mo = i, (54) 

where dy(£) = dRe£ dlm£ is a measure for the integral. The normalization constant 
corresponding to the xy-representation is obtained by an additional factor g/{2^2), 
i.e., 

f(x, y ) = ^=ex.p[~^(x 2 + y2 )]. (55) 

The corresponding normalization is 

Jdxdy /(x,y)/(x,y) = 1. (56) 


Using the creation operator (3' and the relation 

e +«> e -« = e+ «- (z£+I) e -« = (57) 

it is easy to show that the ^-representation of the lowest eneregy state (Lowest 
Landau Level) = (£, £|n = 0, l) is 

//“(«,«) = (58) 

The physical meaning of the quantum number £ is seen by evaluating the expectation 
value of the square of the distance from the origin r 2 = x 2 + y 2 = 8££/ g 2 with respect 
to the lowest Landau levels as 


($$)e 


J d/i(0f e LLL (z, omtf LL & o 

2 e+1 


7r£! 

2^+2 roc 


dm m 




e +1 


dr r 2e+3 e~ 2r2 


(59) 

(60) 
(61) 
(62) 


Here the spherical coordinates £ = r exp(*</>) is used to evaluate the integral. Thus the 
quantum number l corresponds to discrete values of radius for the classical circular 
motion in a constant magnetic field. 

The complete orthonormal basis in the ^-representation (£,£|n,£) can be con¬ 
structed in a similar manner, namely by applying (at) n /Vn\ to / £ iiL (£,£)- It is 
expressed in terms of hypergeometric functions. Such an explicit representation of 
the solution will not be used in the following. 



10 


Jun Suzuki 


§3. Second Quantization of the Dirac Field 


The second quantization for the Dirac held in terms of the complete orthonormal 
basis is achieved as follows. 1 ' 5 ^ 

i’ix) = + dy2(x)e He ” t ), (63) 

A 

where x represents x^ = ( t,x ), and A represents all collective indices, e.g. the 

energy level n, the spin s, and so on. ip A (x) and ip A (x) correspond to the positive 

and negative energy (±e^) solution of the single particle Dirac equation in external 
fields, i.e., 

[7°(±e±) + 7 l pi -e£- m]i/) A (x) = 0. (64) 

Similarly, the conjugate held ^{x) = ijA 7 0 is quantized as 

■$( x ) = 5^(&Ii^i( aj ) e+ * e ” t + d A ijj A {x)e~ l£nt ). (65) 

A 

In the same manner as the held quantization is done in terms of the plane wave 
basis, the equal time anti-commutation relation for the Dirac held is employed to 
satisfy the positivity of the eigenvalues of the Hamiltonian : 

WU x),$ x {x')} =5 pX 5(x-x'), (66) 

X0=x' 0 

where indices p and A denote components of the Dirac spinor. It then follows that 
the annihilation and creation operators b A and 64 for the positive energy particles 
satisfy 

{b A ,b\,} = 8 AjA ', (67) 

and the annihilation and creation operators d A and (A a for the negative energy par¬ 
ticles satisfy 

{d A ,d\,} = 8 a , A ’- (68) 

In order to manifest diagram technics, define the propagator for the Dirac held by 
the vacuum expectation value of chronologically ordered held operators : 

iS(x; x') p \ = (0\T'4> p (x)'$ x (x')\0), (69) 


where T is the time ordering operation. For fermion held operators A and B, the 
time ordering operation is dehned as 


T{A(t)B{t')) = 


+A(t)B(t') (t > t') 
-B(t')A{t) (t < t') 


(70) 


This dehnition can be extend to an arbitrary number of fermion held operators. In 
this case a sign factor (—l) n is assigned with n corresponding to the number of per¬ 
mutations of operators in order to bring the original expression to the chronologically 
ordered expression. It is also convenient to express the time ordering operation as 

T{A(t)B{t')) = &(t - t')A(t)B{t') - e(t' - (71) 
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where 0(x) is the step function. In the above definition of the propagator (69), |0) 
is the vacuum of the Dirac field, i.e., 

MO) = MO) = 0 . (72) 

This vacuum is not to be confused with the vacuum used in the previous section to 
construct the solution of the Dirac equation. The quantized Dirac field ip(x) satisfies 
the Dirac equation : 

(ilj) — m)ip(x) = 0, (73) 

where D. mu = + ieA^. The bare propagator (69) satisfies 

(F0 — m)S(x ; x') = 5(x — x'). (74) 

It should be noted that the propagator depends on two points x and x' in general 
(not the difference x — x 1 ), and hence the Fourier transform of the propagator does 
not take a simple form. 

3.1. Propagator in a constant magnetic field 

Following the above argument, the propagator for electrons in a uniform mag¬ 
netic field is constructed by 

OO « 

VM) = (x,y,z\ E E / dp z (be,n,s(p Z )\^s(i,n,p z -,t)} + 4 ,n,s(P^)l^+ 2 (A«,Pz;*))) • 

£,n=0 s=l,2 d 

(75) 

Here the projection (x, y, z\ is taken in order for the field ip(x) to be defined at 
all space-time points x. As was mentioned earlier, the n = 0 case is allowed only 
for s = 2. Using explicit expressions for the stationary solutions \if a ) (44-47), the 
propagator is 


iS(x;x') = (x,y\^2 C 

e,n J 


ip* C 


2£r, 


0(r)e 


(e n + m)D n (t) C n (£,p z ) 
-C n (£,p z ) - (e n - m)D n (£) 


-\p+ i£nT (^ £n m )Dri(£) C n (£,Pz) 

[ ’ V C n (£,p z ) -(e n + m)D n (£) 


I x',y'), (76) 


where ( = z — z' and r = t — t'. To obtain the above expression, we also used 
the formula (z\p z ) = e lzpz /\/27r and its conjugate (p z \z) = e~ lzPz / y/2n The 2x2 
matrices C n (£,p z ) and D n (£,p z ) are defined by 


/~i ( n \ _ f PzAn—l,n—l(£) i\/r^QA n — l,n(^) 

A ,PZ> ~ \-iyfceK,n-l(Z) PzA n , n (£) 


D n (£) = 


_ fA n —i, n —i(£^ 


0 


0 


(77) 


(78) 


A n , n (£) 

together with a projection operator A nn >(£) = |n, £)(n',£\. The step function 0(x 
has the Fourier transform of the form : 


0(x) 



gikx 

k — ie 


(79) 
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where e is a small positive number which goes to zero at the end of calculations. 
Using this formula, the propagator is expressed by 


l.n 


e ~ipor+ip z Q 
PO ~ £ n + ie 


(PO + m)D r 

-C n 


C n 

-(po - m)D r 


\x ,y , 


(80) 

Using the relation (z\p z ) = e lzpz /V2 tt and its conjugate ( p z \z } = e IZPz /\/2tt, it is 
also possible to write the propagator as an operator in the form of S. Then, the 
usual propagator is considered as a matrix element in space-time coordinates, i.e. 
S(x] x') = (x|Sja/). The operator S is referred to by a Green’s operator to distinguish 
it from propagators in particular representations such as the space-time coordinate, 
or the four momentum representations. The Green’s operator satisfies 


(Jtl - m)S = I, 


(81) 


which is equivalent to eq. (74). This description is more convenient in the following 
discussions. 

It is also convenient to obtain a coordinate representation for the electron prop¬ 
agator, which will be used later. The main concern is to evaluate 

(x,y\A n y(£)\x',y'), (82) 

equivalently, 

OKM')- ( 83 ) 

For this purpose it is easier to carry out the summation over £ first, then to evaluate 
the following function. 

£n,n'(£,<£; £'<£') = £<aK,n'(^r,e>- ( 84 ) 

l 


Since 


(£,£l M) = 




Vinlf! \/2 


the summation over i gives 


2 ,-d + i -d' + Z,' 


ir\/n\n'\ \f2 

2 e-«-^"' 


T(- 


V2 




(-9 + 2e) n (-9' + 2 f) n 'e 2 ^' 


ir\/2 n+n ' n\n'\ 


(85) 


( 86 ) 


(87) 


There are three cases for the values of n and n' : i) n = n' , ii) n — 1 = n! , and iii) 
n + 1 = n!. 

We have for the case i), 

(-<9 + 2£) n {-d' + 2^) n e 2 ^' = {-5 + 2£) n e 2 &\-d' + 2? - 2£) n (88) 

= e 2 &\-8 + 2£ - 2£') n (2£' - 2 f) n (89) 
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3=0 


e m' - 2£') n -i(-d) j {2? - 2£) n (90) 


= 2 n n\e^' ]T 


3=0 


n\ [-2(a-m-?w 


j! 


Use the Laguerre polynomial defined by 


L n(x) = Y[,; 


3=0 


n\ (— x)i 


then 


£»,»«, f; ('(') = - «f) 

7r 

Similar calculation yields for the case ii) : 

£.+!,»(«, (;('(') = - )'> ( 

vr v« +1 v 

where L„ ^ (x) is the associate Laguerre polynomial dehned by 

4 1) w = E^ +lv “ l) ' 

3=0 


(91) 

(92) 

(93) 


^n(2ie-m+4 1 2 1 (2|c-e 


t? 12 


n-jy j! 


(94) 

(95) 


The case iii) is obtained by taking the complex conjugate of eq. (94) together with 
exchanging £ and £ <-* In the same way as before, the following convention 
is understood, 

L_ n (x) = L^ n (x) = 0 for n < 0. (96) 

Also an additional factor in order to get the xy-representation should be remembered. 
Therefore, the explicit coordinates representation of the propagator is 


dp 0 f dp z e -ipoT+ip*C-|£-<fl 2 +fr-£f' 


sW= A U? 


vl - 4 + it 


(p 0 + m)D n 
-C n 


Cn 

-(Po ~ m)D r 


with the definition of 2 x 2 matrices : 


C n = ~ 


Pz-Un—1 


iV2e(^-e)(L n -i+L^_ 2 ) 


*v^e(e-0(^n-i + 4- 2 ) 

Pz U n 


, (97) 


(98) 


= { L V L n ) ■ ("> 

In this expression the argument of the Laguerre polynomials is a function of 2|£ —£ 7 | 2 , 
i.e. L^ 0,1 ^(2|£ — ^| 2 ). Therefore the only term, which depends on two points £ and 
f', appears as an exponential factor exp(££' — ££'). This, of course, agrees with the 
general statement that the electron propagator is gauge dependent object with an 
additional phase factor of the form : exp (ie f , dx^A^fx 11 )). 
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3.2. Photon propagator 

The propagator for a quantized electromagnetic field, the photon propagator, is 
defined by 

-iD^ u (x - x') = (0\T A fl (x)A u (x , )\0), ( 100 ) 


where A^(x) is the second quantized field for photons. Associated with the gauge 
degree of freedom of A /( (x'), the photon propagator is not determined uniquely. The 
most common gauge is the Feynman gauge, which is expressed in Heaviside-Lorentz 
units as, 


-fV ( x -x') = - 


9n v 


4 -7T 2 (x — x') 2 — ie 


( 101 ) 


Its Fourier transform is 




( 102 ) 


In the following, the Coulomb gauge is also employed. The Fourier transform of the 
photon propagator in the Coulomb gauge is 


Doo(k) = -^2 
D 0i {k) = 0, 


Dij(k) = - 


1 


k 2 - k 2 + 


;A°V fc.2 >■ 


ie 


§4. Modified Dirac Equation 


(103) 

(104) 

(105) 


The mass operator J\4 is defined by a summation of all irreducible diagrams 
which contain two external electron lines. Irreducible diagrams are defined by dia¬ 
grams which cannot be reduced to a sum of simpler diagrams by cutting lines. To 
the lowest order correction, it is represented by this diagram 



(106) 


Let S be the bare Green’s operator for electrons, which is constructed directly from 
the solution of the single particle Dirac equation in external fields. Let Q be the 
exact Green’s operator for electrons, then it satisfies the following relation among 
the bare Green’s operator and the mass operator : 




(107) 


In the above diagrams, the thick line represents the exact Green’s operator iQ, and 
each shadowed circle corresponds to the mass operator —iA4. It should be empha¬ 
sized that this is obtained under an assumption that the exact Green’s operator 
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can be decomposed into a summation of irreducible diagrams. Express the above 


diagram relation as an equation by 

G = S + SMS + SMSMS + ••• 

(108) 

= 5 + 5M(5 + 5M5 + ---) 

(109) 

= S + SMG, 

(no) 

thus 

G-S = SMG- 

(HD 


Inserting the following identity, this can be expressed in the coordinate representation 

[dx \x)(x\ = I, (112) 


Q(x;x') — S(x-,x')= dx\ / dx 2 S(x] x\)M(xi\X 2 )G(x 2 ] x'). 


(113) 


Useful relations are derived by applying an operator i/l = f — eA to eq. (Ill) from 
the left : 


UlQ = $S + UISMG (114) 

= (m + I)S+(m + I)SMG (115) 

= mM + I + MG, (116) 

where [l/l — m)S = I is used. Define a new Green’s operator M by 

M=m + M, (117) 

where the coordinate representation is 

M(x; x 1 ) = m5(x — x') + M(x; x'). (118) 

The exact Green’s operator satisfies 

(. Ul-M)G = I • (119) 

In the coordinate representation, 

HpG{x]x') — f dx± M(x] xi)G{xi; x') = 5(x — x'). (120) 


Therefore it is found that the exact Green’s operator satisfies the modified Dirac 
equation (119) or (120), where the mass rn is replaced with the operator A4 = 
m + A!. 16 ) Note that the poles of the propagator correspond to the energy levels of 
the system as seen by the expression (80). Then this modified Dirac equation for the 
exact Green’s operator is particularly important, since the poles of the exact prop¬ 
agator contains information about the exact energy levels of the system. Inverting 
the argument to derive eq. (74), the exact propagator assures the exact solution |^) 
of the modified Dirac equation : 


{Ul-M) |<F> =0, 


(121) 
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or, 

*(x)-/dx'M*;*W) = 0 . ( 122 ) 

Choose the positive energy £ solution 1 Eg), i.e., 

P 0 \Ve)=£\#e), (123) 

then eq. ( 121 ) is rewritten as 

( 7 °£ + fPi - e4 - M)\E £ ) = 0. (124) 

Assume that the exact eigenstate is slightly different from the eigenstate of the 
original Dirac equation due to radiative corrections. Under this assumption the 
eigenstae is written as 

|%) = IlM + \Sip), (125) 

where | V’e) satisfies 

{^°E + 7 l Pj - efi - m)\'ip E ) = 0 . (126) 

By our assumption, AE = £ — E and ||(h/’)l are small. Substituting (125) into 
eq. (124) and using (126), 

(£ - E)\^ e ) = 7 °M|^) - 7°[(7°£ + -e4~ M}\5^). (127) 

Multiply both sides by('0e| from the left, 

£ ~ E = {iPe\7°M\^e) ~ Pi -e£- M]\Sip). (128) 

By egarding the operators in the last term as acting on (i/je\ and using the conjugate 
Dirac equation, 

{iPe\i°[i°{-E) + 7 i P i + ei + m\ = 0, (129) 

eq. (127) becomes 

£ - E = (ipEh 0 M\ipE) + (V , E| 7 °[- 7 °(<f + E) + 2m + 2e^ + M}\8il>). (130) 

The second term in the right hand side is regarded as a next order contribution. 
Therefore the energy shift is found in the perturbation expansion to the lowest order 
as 

£-E = ^ e \i°M\^e). (131) 

In the coordinate representation, this is written 

£ — E= Jd 3 x Jd 3 x' tpE(x)M.(E\x-,x') / ipE(x'), (132) 

where M(E\x] x') is the Fourier transform of the mass operator in the coordinate 
representation with respect to time variable, i.e., 

M(E\x-,x') = f dr e luJT M(x\x') . (133) 

J uj=E 

A time dependence of the form t = t — t' of the mass operator is guaranteed for 
homogeneous external fields case. 
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§5. One Loop Correction to the Mass Operator I 


The one loop correction to the mass operator corresponding to the diagram (106) 
is calculated in this section. To this end it is useful to adopt the mixed representation 
for the propagators, namely momentum space for the time and z-component of 
coordinates, while we use the position coordinates for the x, p-variables. The reason 
for this choice is clear, since the bare electron propagator depends on t — t' and 
z — z!. Therefore, ordinary Feyrnan rules in momentum space is applied for the po 
and p z variables. In order to review notations, a summary for the propagators in 
this representation is given below. The bare electron propagator is 


S(po,Pz\£.,= - V — 2 - 2~i~ — r n(po,Pz\f; ¥), (134) 

K „ Po - £ n + 


where the energy levels are e n = \Jm? + p 2 + g 2 n , and a 4 x 4 matrix T n is 

r,(po,p,\(-i',t-(') = ( (po + ’ 7 )Dn , A (135) 




PzU n — 1 


iV2g(Z - e)(Ln -1 + 4 - 2 ) 

/|2 . _ /L n _ 1 (2ie-ri 2 ) 


iV2 6 (£-i')(L n ^+ 44 ) 


Dn{\t-m = 


-PzLr, 


0 


0 

L n ( 2ie - e 


/| 2 ' 


The photon propagator in Coulomb gauge is 


Doo{ko,k g \Z-?,$-?) = - 


' d 2 k± 


(27r) 2 k 2 

Doi(k 0 ,k z \£, - £',£-£') = 0, 


A #0 Al£-£', <£-<£') = - 


d 2 k ± kik 

"TJT > > 


(136) 

(137) 


(138) 

(139) 

(140) 


(27t) 2 fcg — k 2 + ie v 13 k 2 

where d 2 k± = dk x dk y , k = \/ 2 (k x + ik y )/g, and k = \k\. 

In the hrst approximation, the mass operator is approximated by the diagram 
(106), 


- iM{po,p z \t eo 

= ~ e2 fifr f lur - €',t-t')s(po - ko , pz - k z i^eA'A'b"- 

(141) 


As was shown in the previous section, the energy shift of the ground state due to the 
radiative correction is obtained by the expectation value of the mass operator with 
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respect to the ground state. The ground state (E = m, p z = 0) of the electron in a 
constant magnetic field in this mixed representation is given by 


V’gGCI 2 ) = {Z,£\*l>g) 



o), 


(142) 


where w\ = ( 0 , 1 , 0 , 0 ) is the eigenspinor, and a phase factor exp (—imt) is dropped. 
To reduce the amount of calculations, we also choose l = 0 state, but all other 
choices of the lowest Landau level give the same result. This comes from the fact 
that the ground state degeneracy is not affected by the first order radiative correction. 
Therefore, the energy shift AE g for the ground state is 


AE g 


{^gh°M\^ g }\ po=mpz=0 

I M£) J dug) ^(|£| 2 )^(l£ / l 2 )(?,?l7 0 -M|e , ,£')b= m , Pz =o 


E / Mt) [ e- |?|2 


l«T 





(143) 

(144) 


e -|£-£'l 2 +ft'-£l v / 1 - - n 

(„,_fr,)2_ e 2+ it (t5® 27 r " (m “ k °' ~ kX ~ {, ' ? " ?h ” 2 

+ £ _ (((( w 2 yr n (m - I; : \, c' 7 f I) , (145) 


where e n 


ybr ? 2 + (—k z ) 2 + g 2 n. Changing integration variables to £ and y 


£-£', 


. ^ 4ie 2 
AEg = 2 ~ 

7T 2 


Y f J Mx) 


d A k e~ 2 ^ 2+2x ^~ 2 ^ 2+iKX+iKX 


(27r) 4 (m — ko) 2 — e 2 + ie 


x ( -jj2 W2 ' r n(m - ko, —k z \x, x)tw 2 + 


5ij — kikj/k 2 
kn — k 2 + i 


ie 


wifTnirri - k 0 , -k z \x, x)l 3 w 2 


(146) 


Carry out the integration over £, which is an ordinary gaussian integral, to get 
2 ie 2 f , N r d A k e - 2 lxl 2 +i«x+*«x 




1 


7T 


X] / d Mx) 


(27t) 4 (m — fco ) 2 — £ 2 +ie 


x ( r n {m - k 0 , ~k z \x, x)j w 2 + 


Sij — kikj/k 2 


l’n — k 2 + 


— w 2 i l r n {m- k 0 , —k z \x, x)l Jw 2 


ie 


The matrix elements are evaluated as 

^27 °r n (m - k 0 , -k z \x,X)l°w 2 = - 
Using the formula, 

i fa 0 \ i fa + 2d 

° [o d) a ^ \ 0 


0ko - 2 m)L n ( 2 |x| 2 ). 

° 1 

d + 2aj ’ 


(147) 

(148) 

(149) 
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W 2 'y l U n (m - k 0 , -k z \x,x)l 1 W 2 = -k 0 (L n (2\x\ 2 ) + 2L n -i(2\x\ 2 ))- (150) 

Using another formula, 


kikjcr 1 


i fa 0 


U = k 2 


+ (a - d)k 2 ,kiu\ 


yields 


-y L w 2 'Y l r n (m - ko, -k z \x,x)l J W2 


= --^[fc 2 L n - l( 2| X | 2 ) + k\{L n { 2| X | 2 ) - L n _ 1 (2| X | 2 ))]. (152) 


Hence, 


AE g = 


Y jMx) 


' dJk e - 2 lxl 2 +*«x+*«x i 

(27r) 4 (m — ko ) 2 — e 2 + ie k 2 

ko r ; .2 it , t \ t.2/ 


x '|(fco — 2 m)L n + -j-p —^2 — \k 2 {L n + L n _ i) — k 2 (L n — L n _i)] j> . (153) 
The integration for x is done with the aid of the formula 

r°o —y 2 /4 (n./O' )2n 

/ dx xe~ x L n (x 2 )J 0 (xy) =-- -7—, (154) 


where Jo(x) is the Bessel function of order zero, which has an integral representation 
of ^ 

Jo(x) = — f# e" ia;sin ^. (155) 

27T J-ir 

Going to spherical coordinates x = rexp(fi^), the integral in question is evaluated 
as follows. 

I n = J dyi(x) e ~ 2M2+iKx+i ^L n ( 2|x| 2 ) (156) 

roo 

= 2ir dr re~ 2 ’ 2 L n {2r 2 ) Jo(2|k|t) (157) 

Jo 

= \e- k l/e 2 e n {kl/ e 2 ), (158) 

where fcj_ = «/fc 2 + A; 2 and the function e n (x) is defined by 


Then x integration leads to 


AE g = —ie 




e n {x ) = —r. 

n! 


e -fc 2 7^2 


(27t) 4 (m — /cq) 2 — e 2 + ie k 2 


x < (A: 0 - 2m)e n (k\! q 2 ) + 


fcp [fcje n (fcj_/g 2 ) + (fc 2 + k 2 z )e n -\(k\!g 2 ) 
fcg — k 2 + ie 


. (160) 
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Lastly, the contour integrals to evaluate ko integral are used, 

fdko k 0 - 2 m = m_ , - 

J 27 n (m - ko ) 2 - £ 2 + ie 2e n ’ 

/ dfco fc 0 m 1 

27ri [(m - k 0 ) 2 - £ 2 + ie] [A’q - k 2 + ie] 2e n (e n + A:) 2 - m 2 ’ 

Therefore, the energy shift of the ground state is found in the lowest correction as 


OO « 

AE g = 2irma / 

--n J 


d 3 k e-^/e 2 

(27r) 3 £ n k 2 


k 2 ± \ k ± en + ( k 2 + k 2 )e n -\ 
Q 2 ) + (e n + k) 2 - m 2 


where a = e 2 /(Anhe) is the fine structure constant in Heaviside-Lorentz units, and 
the same convention e_i(fc^/^ 2 ) = 0 should be noted. 

This integral (162) exactly agrees with the one obtained by Luttinger based on 
different arguments. 8 ) Note that this integral is essentially two dimensional since it 
does not depend on the direction of k± = (k x , k y ). Because of its complexity, it is 
unlikely that we can obtain an analytic expression for this integral, especially the 
second term in the square bracket. 

§6. One Loop Correction to the Mass Operator II 

In this section an alternative expression is given for the first order radiative 
correction to the ground state energy. By reexamining Luttinger’s result, it becomes 
clear that his expression is equivalent to Schwinger’s expression. This will be shown 
briefly below. 

6.1. Comparison to the proper time representation 

The photon propagator in the Feynman gauge is employed instead of the Coulomb 
gauge in this section. Since the calculation is almost same as was done for the 
Coulomb gauge, we only show the result. After integrating over holomorphic coor¬ 
dinate variables £ and the energy shift for the ground state is 


AE g = -2ie 2 Y J 


d 4 k _ fc 2 /p2 me n [k\!g 2 ) + k 0 e n -i{k 2 ± /g 2 ) 
(27r) 4 6 [{m - ko) 2 - £ 2 + ie] [k^ - k 2 + ie] 


Using Schwinger’s trick to bring denominators into the exponent, i.e. i/(a + ie) = 
/ 0 °° dae ia ( a+ie ' ) , eq. (163) becomes 

E r d 4 k f°° f°° 

/ 77m / / dQ 2 ( me n + k 0 e n -i) 

n J (2vr) 4 Jo Jo 

x exp{— k]_/g 2 + ia\[(m — ko ) 2 — £%} + i(*2[ko ~ k 2 ]}. (164) 


A small positive e for the integral convergency should be reminded. A summation 
over Landau Levels n = 0,1, 2, • • • can be carried out at this point, which is just an 



QED in a Uniform Magnetic Field 


21 


exponential series; 


AE g = {-if 2e 2 


' d 4 k 


dai / da.2 (to + k^e zaiS ) 


x exp 


02)^0 — 2 iaimko — (1 — e iaie + ia2Q 2 )k]_/g — i{a\ + a.2)k 


(165) 


Lastly, integrals over the four momentum is done to yield 

9 / . - 2 


_ ima 

AE„ =- 

9 2 vr 


J da\ J da .2 


g 2 (ai + a2 + a2e 1Q ) m 2 a? , 

exp (—1 - 

a 1 + a 2 


(aq + «2) 2 (1 — e iai & + ia 26 2 ) 


(166) 


To see the exact correspondence, we change integration variables aq,a 2 to x = 
g 2 a2/2 and u = aq/(ai + a^)- Thus, 




ma f°°dx l + ue 2lx i2m 2 ux. 

— — du- ---z—- -j- exp(-j—) 

27T Jo x J 0 1 — u + ue “sini/i g^ 


(167) 


This is identical to the expression obtained by Schwinger except for the contact term 
in his expression. 12 ^ This contact term will be discussed in Sec. 8. 

It is also straightforward to check eq. (163) is indeed equivalent to eq. (162) by 
integrating over in eq. (163). This leads to 


00 „ 

AE g = 47rma / 

n=o 


d 3 k 


o-kf/e 2 


(27r) 3 (e n + k) 2 — m 2 


(— + \)e n {k\/g 2 ) + —e n -\{k\/g 2 ) 

S’ff K 


which is identical to eq. (162) with the aid of relation e 2 = in 2 + k 2 + ng 2 . 


(168) 


6.2. Alternative representation 

We next derive an alternative representation for the energy shift AE g . The 
derivation is based on use of the following identity in eq. (163), 

AB = I dt [tA + {l-t)Bf (169) 

After rewriting the product of two terms in the denominator of (163) using this trick, 
elementary integrations over ko and k z yields 


AEg 


ma 

47T 2 


00 


E 

n=0 




e n {kj/g 2 ) +te n -i{k\/g 2 ) c _ k 2 /g i 
(1 — t)kj_ + m?t 2 + ng 2 t. 


(170) 


We can further proceed t integral in a simple manner. Define 


I{lo; nrf) 


1 f'd 1 

2 Jo t 2 — {u> — nrj)t + oj ’ 


(171) 


with uj = kj_/m 2 and rj = g 2 /m 2 . Then, depending on the value of u (or kf), t. 
integral /(ca; nrf) can be expressed in terms of inverse of either tangent or hyperbolic 
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tangent as follows. 


I(w; nr)) 


J VDZ Vo > +nr)> 

1 1 tan -1 !-^^-) 

V \T-DZ y u+mi > 


; 0 < uj < ijO- or uj + < uj 
5 UJ — < LO < UJ + . 


Here D' w and ui± are defined as 

Da = lo' 2 — 2 (my + 2)lo + n 2 if = (lo — lo + )(uo — w_), 
w± = (Vn?/ + 1 ± l) 2 . 


(172) 


(173) 

(174) 


Using this definition for I^ n \w,r)), we arrive at the following alternative expression 
for the energy shift : 


. _ ma 

AE„ =- 

9 2 tt — jo 

n=0 J ' J 


OO /»00 




.W. 


2 e n (—) + (w - nij)e n -i(—) 
rj rj 


I(w; nr)) 


+ I 1 °g(!^±l) en _ i (^)[>. ( 175 ) 


Therefore, we have obtained one integral over the transverse momentum \k±\ and 
summation over the Landau levels n. This expression (175) may have an advantage 
compared to the standard double integral expression (167) particularly for numerical 
purpose. 


§7. Anomalous Magnetic Moment 


We would now like to bring our attention to discussion on the anomalous mag¬ 
netic moment of an electron for the weak field case. The magnetic moment g, e of the 
electron is defined by 


dE a 


l^e — 


dB 


(176) 


B =0 


In other words, by applying a weak external magnetic field to the electron, the ground 
state energy is shifted by amount —)i e B. Therefore the main concern is to examine 
the above energy shift in the weak magnetic field limit, or a small g = s/—2eB limit. 
Our derivation of the anomalous magnetic moment is based on the method used 
by Luttinger. 8 ) To appreciate his elegant derivation, we use the expression (162) 
obtained first by him. Eq. (162) is rewritten as 


. _ ma 
AE„ = - 

9 2i r — jo 

n=0 J u 


where 


oo /»oo r 

S/o dk±k B- 


00 dk e- k y^e n (kl/g°A 

— OO 


k 2 


— —I - kj_F n (k±, k z ) + ( k 2 + k z )F n +i(k±, k z ) 

£r). 


Bn(k_ L — 


£n[(£n + k) 2 - m 2 ] ’ 


, (177) 

(178) 
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and A = \Jk\ + k 2 . Note that the function exp(— x)e n (x) (x = k\/g 2 ) is the Poisson 
distribution. When x becomes large, the Poisson distribution is approximated as the 
Gaussian distribution centered at n ~ x with its width ~ 1 /\fx. Then, in the g —> 0 
limit, exp(—x)e n (x) has a sharp peak around n — k 2 L /g 2 . On the other hand, the 
rest (terms inside a square bracket in eq. (177)) in the above integral is a slowly 
varying function of n. Therefore, the above integral is approximated by expanding 
a slowly varying function around the peak. Namely, let f(n) be the slowly varying 
function of n, then the summation is approximated by 


X e n {x)f{n) 

n=0 


™ ™ fU)(x) 

E^ e »wD L ir (B ‘ I) ' 1 

n=0 j =0 J' 

f(x) + xf ( ‘ 2 \x)/2\ + x/® (x)/3! + (3x 2 


(179) 

+ x)/ ( 4 ) (x)/4! H-, 

(180) 


where e n( x ) = exp(x) is used. In order to evaluate the last term in eq. (177), 

the following different version of that summation formula will be used, 


X e n -i(x)f(n) 

n= 1 


e n -i(x) 2_^[ -q—( n-x) J ] 


(181) 


n =1 j =0 J 

f(x) + /^(x) + (x + 1) (x)/2! + (4x + l)/®(x)/3! 

+ (3x 2 + llx + 1)/ (4) (x)/4! + • • • . (182) 


Since the magnetic moment is given by the term proportional to g 2 = —2eB, it is 
enough to keep only those terms in the following. Using this resummation technique, 
the terms proportional to g 2 are obtained by 


^ _ x , . 1 3f? 2 

> e e n {x) -> ■ 

n =0 £n 

oo 

5 ~2e~ x e n (x)k 2 ± F n 

n =0 

P 2 k 4 : 


k 2 


8 (m 2 + A : 2 ) 5 / 2 ’ 


+ 


16 A 2 (m 2 + A 2 ) 1 / 2 \ {m 2 + k 2 ) 2 A 2 (m 2 + A 2 ) 


and 


(183) 

(184) 


k 2 [{m 2 + A 2 ) 1 / 2 + A :] 2 


X e n -i(x)(k 2 + k 2 z )F n 

71=1 

Q 2 Aj_(A 2 + A 2 ) f 
16 k 2 (m 2 + A 2 ) 1 / 2 \ 


g 2 k 2 + A 2 
4 k 2 (m 2 + A 2 ) 3 / 2 

3 1 

{mi 2 + A 2 ) 2 k 2 {m 2 + A 2 ) 


(185) 


1 

A 2 [(m 2 + A 2 ) 1 / 2 + A] 



After these resummation, the integral yields the terms proportional to g 2 as 


AE^ = 


mag 

2vr 


dk± 


r -4 


3A 2 


A 2 + k 2 


4(m . 2 + A 2 ) 5 / 2 4A 2 (m 2 + A 2 ) 3 / 2 
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k ± 

8k 2 (m 2 + k 2 ) 3 / 2 


__ 

8 k 2 (m 2 + A: 2 ) 1 / 2 [(?n 2 + A : 2 ) 1 / 2 + A :] 2 


Lastly the following elementary integrals are carried out : 


(186) 


and 



kj_ 

4 

(187) 

k 2 (m 2 + A 2 ) 5 / 2 

9 m 2 ’ 

k±(k 2 + k 2 ) 

8 

(188) 

k A (m 2 + A 2 ) 3 / 2 

3 m 2 ’ 

kj_ 

4 

(189) 

fc 4 (m 2 + A 2 ) 3 / 2 

3 m 2 ’ 


dkl r ik _ 4 = j_ r ix 1 = j_ 

J oo ~ k 4 (m 2 + A 2 ) 1 / 2 [(m 2 + A 2 ) 1 / 2 + A] 2 3m 2 7 0 (1 + a:) 2 3 m 2 ’ 

(190) 

where the integration variable is changed to k/(m 2 + k 2 ) 1//2 = x in the last integral. 
Combining these integrals, the final result for the energy shift in the weak magnetic 
field limit is 



AE^ 


ag 2 

1 CO 

i 4^ 

2irm 

\4 9 

a 

f-eB 

~2n 

\ 2m 

a 

~2tt' 

HbB, 


18 14 

4 3 + 8 '3 



(191) 

(192) 

(193) 


where ^b = |e|fo/2mc is the Bohr magneton, and all other terms except those pro¬ 
portional to q 2 are neglected. The geomagnetic factor g of the electron is 2 from 
the Dirac equation. If the true geomagnetic factor including radiative corrections is 
expressed as g = 2(1 + a), then the interaction energy between the magnetic mo¬ 
ment and external magnetic fields B is E = —(1 + o)^bB. Hence the energy shift 
is written as AE = —agsB. From the above result, a is found in the first radiative 
correction as 


_ 9 ~ 2 _ a 
2 + 2t r' 


(194) 


This is precisely the result first obtained by Schwinger with the aid of subtractions 
to remove the divergence of integrals. 17 !’ 18 ) 

From the above derivation, it now becomes clear that the physical meaning of 
the anomalous magnetic moment and the role of Landau levels to it. For the weak 
external magnetic field case, separation of these Landau levels are sufficiently small 
such that electrons can be excited within a very short time and come back to the 
ground state by interacting with photons surrounding them. As we have seen, in 
quantum electrodynamics these virtual transitions happen equally even to the higher 
Landau levels. This is understood from the fact that the Poisson distribution is well 
approximated with the Gaussian distribution in the weak magnetic field limit. In 
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this sense, existence of the anomalous magnetic moment is a proof of underlying 
infinite number of Landau levels, even if we do not observe them directly. 

§8. Divergent Terms and Zero Magnetic Field Limit 


We now examine divergent terms appeared in the weak magnetic field limit of 
the ground state energy. These terms are independent of this applied external field, 
i.e., terms of order g°. Using the resummation technics in the previous section, we 
obtain 

. _r„o\ rna f°° „ 

AEt U = -2 vi ik - 

(195) 

mot 
~ ~2 tt 

Obviously, the second term in eq. (196) contains ultraviolet divergent term. This 
divergent term has the form 3ma/(47r) hr(A; 2 iax /m 2 ) (fc max —> oo), which coincides 
with the same divergent behavior as the case of radiative correction to the electron 
self energy without external fields. 1 ) We remark that this divergent behavior is also 
clearly seen in the last term of the expression (175). 

Therefore, the mass of the electron needs to be renormalized. In other words 
the bare mass m appeared in this formalism is not a physical mass. Alternatively, 
the mass should be defined in the limit B —> 0 in the modified Dirac equation (121) 
in Schwinger’s language, since the physical mass is not m but A4. 

We next evaluate terms proportional to g 4 , i.e. next order in the weak magnetic 
field limit expansion. Using the formulae (180,182), we obtain 


\Jm 2 + k 2 ^ 3 1 ( \Jm 2 + k 2 + k 

\Jm 2 + k 2 + k 2 l x/m 2 + k 2 — k 


(196) 


fc^O 


/: 


dk z 


k _l 
k 2 


+ 


Lx/ ml 2 + k 2 \Jm 2 + k 2 iy\Jm 2 + k 2 + k)\ 


OO 

^2e~ x e n (x)fi(n) -> x/ (3) (x)/3! + 3x 2 / (4) (x)/4! (197) 

n =0 
oo 

S ^e~ x e n -i{x)f 2 {n) -> / (2) (x)/2! + 4x/ (3) (x)/3! + 3x 2 / (4) (x)/4!, (198) 

n —1 


where /i(n) is either l/e n or k“j_F n (k±,k z ), and f 2 (n) = (k 2 + k^)F n (k±, k z 
ing k = \/k 2 + m 2 , straightforward calculations lead to 


d 3 ( i \ 


15£ 6 


dn 3 

d 4 
dn 4 



n=kf/g 2 


n=k^/e 2 


8 F ’ 
105^» 8 
16 k 9 ’ 


Denot- 


(199) 

( 200 ) 


and 


d 2 F n 


dn 2 


n=kf/g 2 


3 g 4 _J f _ g 4 

8 k 2 k 5 8 k 4 k 3 8k 4 k(k + k) 2 ’ 


(201) 
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d 3 F n 

On 3 

d 3 F n 

dn 4 


15 g 6 3g 6 

n=k 2^/ g 2 16 k 2 k 7 16 k^k 5 ' 

105^ 8 f 45^ 8 ( 9g 8 

n=fc 2 /e 2 ~ 32k 2 k 9 + 32 k A k 7 + 32 k 6 k 5 

3g 8 3 g 8 3g 8 

32 k 8 k 3 + 32k 8 k(k + k) 2 + 16k 7 k(k + A;) 3 ' 


( 202 ) 


(203) 


Careful evaluations for the above integrals show that there are also divergent terms 
appeared. To identify leading behavior of divergence, we keep finite boundaries for 
integrals. The final expression is then, 


where we have defined the integrals I® (i = 1 ~ 5) as follows. 



AE (e 4 ) = r ^l. 

9 2ir 


— !) _ —l ( 2 ) + —l( 3) _ 293 /(4) + 863 j(5) 

64 672 420 16800 16800 


jd) 
/( 2 ) 
j(3) 
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fc 4 

fc 2 
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1 


+ 


+ 


k 2 k 3 k 2 k(k + k ) 2 &&(A; + A;) 3 



1 

fcfe(A: + fc) 3 ‘ 


(205) 

(206) 

(207) 

(208) 

(209) 


It is easy to see to observe that all integrals except /d) have finite values in the limit 
fc min —i► 0 and fc max —> oo. Those values are 


j(!) = 2 j(2) = 2 j(3) = j(4) = _ 1 


35m 4 ’ 


15m. 4 ’ 


3m 4 ’ 


m 


4’ 


and I 1 ' 5 ' 1 i 


is 


1 

2\Jk 2 + m 2 i 

f k \ 

m 4 

k + Vk 2 + m 2 

V A; + Vk 2 + m 2 J 


J k„ 


( 210 ) 

( 211 ) 


Thus, infinity arises in the limit k m \ n —* 0 as ln(fc m i n /m). However, we notice that 
this infinity is originated from improper use of the resummation formulae (180,182). 
Indeed, these formulae fail for very small k, where the function f(n) is no longer 
slowly varying function of n. The same difficulty is also encountered when one 
tries to evaluate the weak magnetic field limit using the proper time representation 
eq. (167). 12 ^ However, we remark that the anomalous magnetic moment of an elec¬ 
tron can be obtained in the weak field limit without any divergence by employing 
the definition (176). 
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§9. Lowest Landau Level and Strong Magnetic Field Limit 

In this section we utilize the lowest Landau level approximation in the strong 
magnetic field limit. For sufficient large values of B, energy levels are separated 
by large enough such that higher Landau levels will not contribute to the radiative 
correction compared to the lowest Landau level. In this approximation we obtain 
the leading term from the expression (175) as 

AE hhh = ma r 00 ^ 2e -u/ VI ^. ( 212 ) 

2vr Jo 

It is convenient to change integration variables to get 

AEg Lh = !^. e ~Vv (J*dz z e -( 2cos *)A? + J™dz z e -( 2 «>sM/^ . ( 2 13) 

This integral representation shows that, in this approximation, there is no singular 
behavior in the strong magnetic field limit. Further corrections can be obtained by 
taking into account contributions from the higher Landau levels. 

§10. Conclusion and Discussion 

We have shown that a proper choice of the Fock space would not lead to divergent 
integrals when obtaining the radiative correction to the anomalous magnetic moment 
of the electron. However, m and e appeared in this formalism are not physical 
observed mass and charge respectively. These quantities should be defined self- 
consistently. 

Another issue which we have not discussed in details concerns applications to 
strong external field problems in quantum electrodynamics. There have been many 
activities on this issue in the past. Essentially the same quantization scheme is used 
to deal with intense external fields problems. From our discussions the reason is 
clear why we need to quantize electron fields in classical external fields. If one starts 
from a Fock vacuum for non-interacting electrons, as is done in ordinary quantum 
electrodynamics, then one cannot obtain the true vacuum for interacting electrons 
in strong fields in a general perturbative manners. Therefore it is crucial to choose 
the initial Fock vacuum as close as possible the real interacting electrons in external 
fields. 

In this paper we have also derived the alternative expression for the energy shift 
of the ground state in a uniform magnetic field. Detailed analysis of our expression 
will be studied together with applications to strong magnetic field separately. 19 ) 
These are especially relevant to recent astrophysical observations. 

In conclusion, the renormalization procedure is necessary in current quantum 
field theories, since we do not know how to solve interacting systems in general. The 
structures of quantum field theories themselves, however, seem to have less to do 
with the reason for divergent numbers. Rather, the problem is originated from the 
simplifications in our theoretical models, such as an assumption of non-interacting 
particles, infinite volumes, infinite number of particles, and so on. We then have a 
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ground to conclude that a better choice of the Fock vacua may not necessarily result 
in divergence difficulties. 
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